Effect of the spin-orbit interaction 
on the band structure and conductance of quasi-one-dimensional systems 
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We discuss the effect of the spin-orbit interaction on the band structure, wave functions and 
low temperature conductance of long quasi-one-dimensional electron systems patterned in two- 
dimensional electron gases (2DEG). Our model for these systems consists of a linear (Rashba) 
potential confinement in the direction perpendicular to the 2DEG and a parabolic confinement 
transverse to the 2DEG. We find that these two terms can significantly affect the band structure in- 
troducing a wave vector dependence to subband energies, producing additional subband minima and 
inducing anticrossings between subbands. We discuss the origin of these effects in the symmetries 
of the subband wave functions. 
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I. INTRODUCTION 



In 1986, the first experimental realisation of a quasi- 
one-dimensional electron system (QIDES) in a dynami- 
cally confined two-dimensional electron gas (2DEG) was 
achievedtl. Since then there has been an extensive theo- 
retical and experimental effort into understanding their 
basic properties (see, e.g., reviewB and bookcl and ref- 
erences therein). The interest in these systems stems 
from two facts. First, the effective transverse size of a 
QIDES can be easily controlled and made remarkably 
small, down to the de Broglie wavelength of an electron. 
This makes it possible to realise experimental systems 
which have an arbitrary number of occupied transverse 
modes. Second, the high purity of 2DEGs grown by 
molecular beam epitaxy enables the almost coUisionless 
motion of an electron through an experimental QIDES. 
The coexistence of these two factors has made QIDES 
unique objects for the investigation of transport phenom- 
ena yielding, in particular, tke-observation of the ballistic 
quantisation of condwjCtanceQ'El and the so called 0.7 con- 
ductance structureOil. 

The process of electron transmission through a QIDES 
involves the redistribution of incoming electron flux 
among its discrete eigenstates followed by adiabatic 
transport through them. Therefore the determination 
of the electron eigenstates of a QIDES is an integral and 
very important part of the more general quantum trans- 
port problem. This statement is especially relevant to the 
ballistic transport regime where the total conductance of 
a system is completely defined by the number of prop- 
agating electron modes which in turn can be uniquely 
calculated from the energy spectrum and the Fermi en- 
ergyBBa. 

Clearly, the energy spectrum of electrons crucially de- 
pends on the effective geometry of a QIDES as well as on 
external and internal fields acting on them. Among the 
possible internal forces, one of the least understood ex- 



amples is the interaction between orbital and spin degrees 
of freedom of an electron: the spin-orbit interaction, also 
referred to as the spin-orbit coupling. Although this in- 
teracjtipp-has an essentially relativistic nature (see, e.g., 
Refs.tHll3), it nevertheless can give rise to an jObaervable 
modification of semiconductor band structureEJcil. 

A quite general theoretical approach to the description 
of tbe||-S|0 interaction is to use the following Hamilto- 
nianll30 which stems directly from the quadratic in v/c 
expansion of the Dirac equation: 
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Here Mq is the free electron mass, p is the momentum 
operator, e is the elementary charge, & = {ax,cry,az} 
is the vector of the Pauli matrices, V{r) is the poten- 
tial energy of the particle, and V stands for the spa- 
tial gradient. The convenience (or universality) of the 
Hamiltonian (|l|) is that it does not restrict one to any 
particular form (model) of the potential V{r), but allows 
freedom of choice depending on the nature and the sym- 
metry of forces present in a given medium. Its form is 
such that it can remove the spin degeneracy in electron 
band structure whilst not actually producing an overall 
spin polarization. 

In a bulk (3D) crystalline environment, the energy 
V{r) arises exclusively from the periodic (microscopic) 
crystal potential. Most multicomponent semiconduc- 
tors have either zincblende (GaAs and most III-V com- 
pounds) or wurtzite (II- VI compounds) lattice structure, 
both of which lack inversion symmetry. DresselhausE2l 
has shown that this property eventually leads to a SO- 
induced splitting of the conduction band into two sub- 
bands. The magnitude of the splitting is proportional to 
the cube of the electron wave number k. 

In MOSFETs and heterostructures the host crystals 
cannot be treated as ideal 3D systems, because the crys- 
tal symmetry is broken at the device interface where 



2D electron or hole gases are dynamically confined in 
a quantum well. The reduction of the effective dimen- 
sionality lowers the symmetry of the underlying crystals 
and results in an additional (linear in k) term in the 
Dresselhaus splitting. Moreover, if the quantum well is 
suffj cjentlt yj-aarrow, then the linear contribution is dom- 
inanlll3'E2rE3 and, e.g., may reach ~ 0.3 ~r-p.4 meV in 
180-A-thick modulation-doped GaAs wellal3. Theoret- 
ical argumentsE^ suggest that this can also be true in 
strained III-V crystals and hexagonal II- VI compounds. 

Along with the microscopic crystal forces, there is an- 
other source of the potential energy T^(r) in 2D sys- 
tems. It is caused by the interface electtic, field that 
accompanies the quantum-well asymmetryaEZI and is di- 
rected along the normal to the device plane. Since 
the typical well width is ^^ 1-10 nm, the interface 
potential turns out to be slowly varying on the scale 
of the lattice parameter and can be considered macro- 
scopic, as opposed to the rapidly oscillating atomic 
field. The mechanism of the SO interaction originat- 
ing from the interface field was first introduced by 
Rashba in RefBa. It also manifests itself as a linear 
in k splittingc^ of the 2D bandstmcture. In a liariety 
of systems induing Si-M0SFET!£3, InAs/GaSbB and 
AlSb/InAs/Al^£3 quantum wells, InGajAs/InAlAs het- 
erostructuresE3Lil, GaAs electron gasesEil it can be made 
to dominate the Dresselhaus terms indicating the signif- 
icance of macroscopic potentials in producing observable 
SO effects in low-dimensional systems. 

In modern nanotechnology there exists a number of 
methods for creating QIDES from 2DEGs: the split- 
gate technioue; wet and dry etching; and cleaved edge 
overgrowthEn In essence, they all exploit the confine- 
ment of the lateral (in-plane) motion of electrons (holes) 
by some transverse potential. Any such potential must 
be essentially non-uniform in space in order to force the 
charged particles to remain within a confined area. The 
spatial variation scale of the confining potential crucially 
depends on the particular fabrication method and varies 
over a wide range: '--^ 10 - 1000 nm. Thus, in QIDES 
one finds another example of macroscopic potentials, viz 
the lateral confining potential, which is absent in higher- 
dimensional structures. The spatial non- uniformity of 
the confining potential gives rise to an additional (in- 
plane) electric field in the system. If the confinement 
is sufficiently strong (narrow and deep), then this field 
may not be negligible in comparison with the interface- 
induced (Rashba) field. Moreover, in nearly square (i.e. 
symmetric) quantum .afdls where the Rashba field is es- 
sentially suppressec£30, the in-plane ( "confining" ) elec- 
tric field is likely to be dominant. This suggests the pos- 
sible importance of the lateral confining potential to the 
SO Hamiltonian (Q) in QIDES. We are not aware of any 
experimental evidence or measurement of the strength 
of the SO coupling resulting from such a confining po- 
tential and cannot therefore quote a grounded estimate 
for the corresponding energy modification. Nevertheless, 
the above arguments seem sufficiently strong to indicate a 



possible new mechanism for the SO interaction in QIDES 
and point towards new transport effects. The existence 
of an additional (and easily controllable) source of the 
SO coupling could catalyse experiments on quantum-wire 
based devices that exploit boihpthe charge and spin of an 
electron, e.g. spin transistorsoEj and active spin polaris- 
ersEj. 

Our belief in the importance of the SO interaction in 
semiconductor QIDES is strongly supported by the fact 
that observable manifestations of a SO-induced energy 
splitting have already been found experimentallyEj in 
another type of QIDES, viz in electron gases trapped 
by dislocational defects in silicon crystals. The source 
of the SO coupling in these systems is an electric field 
perpendicular to dislocations and the SO-related energy 
splitting is linear in the wave number fc, which makes a 
physical analogy between the electron gas on dislocations 
and semiconductor QIDES quite close. The SO interac- 
tion manifests itself in spin-dependent electron conduc- 
tivity along dislocationjLJvhich can be ascribed to com- 
bined spin resonanceESEj of electrons corresponding to 
transitions between SO-split energy levels. 

In this paper we analyse theoretically the effect of the 
spin-orbit interaction on the energy spectrum and con- 
ductance of a long QIDES. As a model for the QIDES, 
we use a strictly 2D electron gas subject to a trans- 
verse electrostatic confining potential. To decide on a 
reasonable shape of the confining potential, we assume 
a sufficiently small effective width (< 300 nm) for the 
QIDES and a low electron concentration (< 10^^ cm~^). 
Combined together, these two factors prevent the con- 
fining electric field from being significantly screened by 
the electron gas. Under these conditions the confining 
potentiaLcap-be accurately approximated by a parabolic 
potentialE3c2l. This conclusion is very favourable to our 
problem because exact analytical expressions for the en- 
ergy spectrum and wave functions of a 2D electron gas 
in a parabolic potential are well known and provide us 
with a good zero approximation in dealing with the SO 
coupling. 

We include the SO interaction via the Hamiltonian (|^). 
We assume that the potential T^(r) which is responsible 
for the SO coupling consists of two contributions: (i) a 
parabolic confining potential with a gradient (or the ac- 
companying electric field) which lies in the plane of the 
2DEG; (ii) a potential which arises from the asymme- 
try of the quantum well with the corresponding electric 
field (Rashba field) being uniform and directed perpen- 
dicular to the device plane. We neglect the crystal-field 
(Dresselhaus) contribution to V^(r). 

The goal of this paper is to reveal the qualitative hall- 
marks of the SO interaction in QIDES rather than to 
construct a complete and realistic spectral and trans- 
port theory. Therefore we use two major simplifications. 
First, we neglect the Coulomb interaction between elec- 
trons. At first sight one may not expect this approxima- 
tion to work in low-concentration 2D electron systems 
where the Coulomb energy may exceed the kinetic en- 



ergy by an order of magnitpde. However, it has recently 
been shown theoreticallycil that the effect of electron- 
electron interactions on SO coupling in such systems can 
be plausibly taken into account via a renormalisation of 
the SO coupling constant. More specifically, this renor- 
malisation leads to an enhancement (by 10 - 50 %) of 
the strength of the SO interaction, which emphasises the 
significance of the SO-related effects in low-dimensional 
electron systems. Our second simplification is the_ ex- 
ploitation of one-band effective mass approximatioro for 
the Schrodinger equation. Within this approximation, 
the infiuence of the crystal forces on electron dynamics 
in the conduction band is reduced to the renormalisation 
of the electron mass and all interband transitions are left 
out.r-.pcspite the obvious simplicity, this approach works 
wellocS for a wide range of semiconductor materials. 

Section || is the central part of the paper. It is de- 
voted to the solution of the problem that we have out- 
lined above a nd to the analysis of the results obtained. 
In subsection II A we introduce the Hamiltonian of a 
2DEG that includes the parabolic confining potential and 
the SO-intcraction term p). We take into account only 
those contributions to the SO Hamiltonian that arise 
from macroscopic (i.e. relatively smooth) potentials, viz 
from the quantum-well (Rashba) and the parabolic con- 
fining potentials. Since the operator (nf) contains the 
Pauli matrices, we seek electron wave functions in the 
form of two-component spinors. Within such a represen- 
tation the Schrodinger equation turns out to be a system 
of two coupled differential equations with respect to the 
spinor components. In next two subsections we solve this 
syst em f or two physically different situations. Subsec- 
tion II B deals with the case of zero Rashba SO coupling, 



where the entire SO interaction comes from the parabolic 
confining potential. Here the Schrodinger equations de- 
couple and the electron wave functions and energy spec- 
trum a re fo und in an explicitly analytical form. In sub- 
section [I C we consider a more general situation where 
both mechanisms of the SO-coupling are present. In this 
case we calculate the electron wave functions and en ergy 
spectrum numerically using the results of subsection II B 
as a convenient basis. Afterwards we analyse the energy 
spectrum obtained in order to reveal the basic features 
of both SO-coupling mechanisms. In subsection II D we 



discuss possible manifestations of the SO interaction in 
the dependence of the ballistic conductance of a QIDES 
on the Fermi energy. Section III summarises the results 
of our research. 



II. THEORY AND ANALYSIS 

A. Problem statement. Spinor equation 

Within one-band effective mass approximationcS the 
Hamiltonian of a QIDES can be written as 



H = 



11 

2M 



VLc(r) + Hso, 



(2) 



where r = {x, y} is a 2D position vector and JLf_is the 
effective electron mass. In line with the Refs.LflLj, the 
lateral confining potential V^ci''^) is approximated by a 
parabola 



VLcir) = 



Muj^ 



(3) 



The quantity uj controls the strength (curvature) of the 
confining potential. The in-plane electric field 'EiLci^) ^-s- 
sociated with Vlc{t^) is given by ELC'(r) = — VVLc(r) = 
-Mw^x. 

We assume that the SO interaction Hamiltonian Hso 
f|I|) in Eq. (pi) is formed by two contributions: Hso — 



TJOt 

J^so 



H 
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The first one, 



^SO' 



arises fronr the asym- 
metriLjoLthe quantum well, i.e. from the Rashba mecha- 
nismEaa of the SO coupling. Since the interface-induced 
(Rashba) electric field can reasonably be assumed uni- 
form and directed along the z-axis, the term H^q can be 
described by the following expression: 



H^o - -^ (^ X P 
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(4) 



The SO-coupling constant a takes values wi th i n .. a f tang e 
(1 - 10)xlO-i" eVxcm for diflFerent systemsBElO-BEl. 
For brevity, in what follows we will refer to the Rashba 
mechanism of the SO coupling as a-coupling. 

The second contribution H^ to Hso comes from the 
parabolic confining potential (H): 



H 



so — 



/3 X 
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(5) 



Here l^ = (h/Mujy^'^ is the typical spatial scale asso- 
ciated with the potential Vlc @)- In Eq. (||) we have 
introduced the SO-coupling constant /3. Comparison of 
typical electric fields originated from the quantum-well 
and lateral confining potentials allows one to conclude 
that a plausible estimate for (3 should be roughly 10% of 
a. Moreover, in square quantum, wells where the value of 
a is considerably diminishedEj'U (by an order of magni- 
tude) the constant /3 may well compete with a. Hence- 
forth we adopt the term /3-coupling for the mechanism 
of the SO-interaction arising from the lateral confining 
potential (y). 

Our objective is to find eigenvalues and eigenfunc- 
tions of the Schrodinger equation H^ = E'^ with the 
Hamiltonian H given by Eqs. (|2|) - (H). The wave func- 
tion ^I* — 'l'(r) = {^|(r) \I'|(r)} is a two-component 
spinor and the energy E is measured from the conduc- 
tion band edge. It is easy to see that the Hamiltonian 
(0) - (P) is translationally invariant in the j/-direction. 
We therefore seek solutions to Schrodinger's equation in 
the form of plane waves propagating along the y-axis, 
i.e. *ix(r) = exp{ikyy)^^i{t). Here t — x/l^ is the 



dimensionless transverse coordinate, ^ix(*) i^ ^^^ trans- 
verse wave function, and ky is the longitudinal wave 
number. After substituting this representation into the 
Schrodinger equation, we obtain a system of two differ- 
ential equations with respect to the spinor components 






(6) 



_ exp[-{tTtp/2)y2] 
\/2"n!7ri/2 



= 0,1,2,... 



(10) 



Here Hn{t) is the Hermite polynomial of n-th order. The 
real wave functions (t)^{t) and (/>?(0 form complete sets 
with respect to the discrete quantum number n and are 
normalised by the following conditions: 









(7) 



(8) 



{kxluY is the dimensionless transverse en- 



where Sx 

ergy, k1 = {2M/'h^)E — ky, and the prime denotes a 
derivative with respect to t. The lengths la and Ip de- 
fined by 



la = h^/2Ma, 1(3 = n^/2MP 



(9) 



are characteristic spatial scales associated with the a- 
and /3-couplings, respectively. We note that the functions 
$11 (i) depend on three dimensionless external parame- 
ters: lui/la, L/lp, and kyl^. 

The equations (g) and (H) are arranged in such a way 
that all the terms which couple them together are col- 
lected on the rhs'. It is interesting that this arrange- 
ment separates explicitly the a and P mechanisms of the 
SO interaction. Indeed, the a-terms enter only the rhs', 
while all the /3-terms are contained on the Ihs'. This 
suggests that the /3-coupling is responsible for forming 
independent ("unperturbed" or "non- interacting" ) wave 
functions ^^(t), while the a-coupling mixes them to- 
ether to form the solution of the whole system (||) and 



B. Zero a-coupling (square quantum well) 

We start the analysis of Eqs. (||) and (^ with the 
situation where the quantum well is square so that the 
interface-induced electric field vanishes. In terms of the 
SO-interaction this means that the Rashba mechanism 
(i.e. the a-mechanism) of the SO-coupling can be omit- 
ted. In other words, the characteristic length la of the 
a-interaction tends to infinity: la —* oo, or l^/la —^ 0. 
The study of this relatively simple case will provide us 
with a convenient basis for treating the general problem 
with finite la- 

Once luj/la has been put equal to zero, the rhs' in 
Eqs. (0) and (R) vanish andrtiiey transform into two inde- 
pendent Hermite equationscJ whose exact solutions are 



fJ0';'^)= / di0f^(i)0^x(i) = <5^„, (11) 



where both ipT] and (j)"^, are taken at a same value of ky. 




FIG. 1. The spinor components (/f>t(i) and (/>j.(t) for n = 0, 
L/^/3 = 0.1, L/L = 0.0, kyL = 10.0. 



From Eq. (HG) one can see that the both "up" {(f)^) 
and "down" (0?) spinor components have exactly the 
same shape as functions of t but are spatially displaced 
with respect to each other by an amount of tfs (see 
Fig. ri). According to the definition (0), this displace- 
ment is a direct consequence of the /3-coupling because 
it does not vanish as long as li^/l/s ^ (except where 
kyluj = 0). At finite values of tp, the complete sets of 
functions <j)^{t) and (f>^{t) turn out not to be mutually 

orthogonal, i.e. {'PTil't'^lr) 7^ f^mn- Instead, the scalar 



cross-product ((t>t\\4'7t ) in the case of sufficiently weak 
/3-coupling {tf} <^ 1) is governed by the following asymp- 



totics: 



{^Wu)-^-oitl), 



(12) 
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±j!i±lt, = ±^^^(MJ, 
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'/'ui-^rr ') ^ T\/^ t, = T\i^ '-fskyi^), (14) 
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0{tPg) for p>2. 



(15) 



The "displacing effect" of the /3-coupUng on spinor wave 
functions [see Eq. HG and Fig. ^] is superficially similar to 
the effect of a perpendicular magnetic field on a QIDES. 
However, the essential difference is that a magnetic field 
shifts both spinor components (/)?(i) and 4>^(t) as a whole 
by an amount proportioM|l to the strength of the mag- 
netic field (see, e.g., Ref.llil), but it does not affect their 
mutual spatial distribution. In contrast to this, as can be 
seen in Fig. m, the /3-coupling causes a spatial separation 
of the spin-polarised electron states 0"(i) and 0"(O- 



a-coupling and therefore to finite values of the coupling 
constant lu/la- For arbitrary non-zero values of l^j/la, 
eigenvalues of Eqs. (g) and (0) cannot be obtained as 
simply as they were in the case of li^/la — 0- To anal- 
yse the energy spectrum in the presence ofj-a-coupling, 
we will follow the scheme proposed in Ref.c3. Namely, 
we first expand the unknown functions $|x(i) in terms 
of the unperturbed solutions 4''\i{t) ( [lO| ) of Eqs. (||) and 

(§■■ 
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(18) 




FIG. 2. The transverse energy ei vs. kyl^j for 1,^/lfj — 0.1 



Although the /3-coupling shifts the spinor components 
(/)?, (i) apart, it nevertheless is not capable of lifting their 
energy degeneracy. The dimensionless transverse energy 
Ex corresponding to the both eigenfunctions (p't, (t) (|l 
is given by (see Fig. 0) 
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(n, ky) = 2n+l- {tp/2y 



(16) 



The total electron energy E then forms parabolic sub- 
bands for each n-th transverse mode: 



E = En{ky 



^ ^ (0) 



(n, ky) 
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For zero /3-couphng {luj/lp = 0) formulaS|-|( 
describe the well-known electric subbandaj' 



(17) 
and M) 



<fuW = E5uCW- 



(19) 



m=0 



We then substitute the expansions (y_8|) and (19) into 
EqSj-dq) and (M) respectively and make use of the prop- 
erty 



dt 



4.- it) 



^n'it) 



b-t\t) 



(20) 



to remove the derivative with respect to t on the rhs' 
of Eqs. (||) and (|^). We next multiply the equations 
obtained by (fi'Ht) and 4>^{t) respectively and integrate 
them over t from — oo to oo. Owing to the orthogonality 
condition (|ll|) , the summation over m is removed by the 
delta-function (5„„i and we find the simple vector rela- 
tions: 



f| = [/| f;. 



gi = t^i 



(21) 



(22) 



Here f|| = {f^^, f\^, . . . J\\, . . .} and gi^ = 
{5? 1 7 St I ; ■ • ■ J 9^\ I ■ • ■ } are vectorised coefficients of the ex- 
pansions (|l^) and (|l9|) respectively, and C/|x are tridiag- 
onal matrices defined by their elements: 
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k 1 

rxjyii^ 



j-j-n.n+l I 'w I n + 1 



e^°\ny 



1/2 



£x- £x yn) 



(23) 



(24) 



C. Finite a-coupling (triangular quantum well) 

We now examine the situation where the effective ge- 
ometry of the quantum well is such that the interface- 
induced electric field is non-zero. An example is the 
triaHjgHlar well which forms at a hetero junction inter- 
faceSO. This interface-induced field gives rise to a finite 



U. 



n 



L f n + l 



1/2 



ex~e^x\n+l)' 



(25) 



Note that the matrices U]^ are neither symmetric nor 
anti-symmetric. 

Eqs. (pi|) - (p5| ) illustrate the effectiveness of the repre- 
sentations ( p^ and (19). Indeed, by using the functions 



(^11 (t) as expansion bases, we have reduced diflFerential 
operators on the Ihs' of Eqs. (ph and (M) to the scalar 
factor Ex — Sx (n). In other words, matrices which were 
supposed to act on f| and g| in Eqs. (|2|) and (E^ turn 
out to be diagonal within the representations ( [iSf ) and 

Our next step is to find the relationship between the 
vectors f|x and g||. To do this, we equate the rhs of 
Eq. (|l^) to that of Eq. ( pj|) and take the scalar product 
of the resulting equation with <j>V'{t). After making use 
of the orthogonalisation condition (O) , we find that 



where the matrix W| is defined by 



(^uI€t)' 



wrr 



(26) 



(27) 



According to Eq. ( |l0| ) [see also Eqs. (12) - (p^j the ma- 
trices W^i are not diagonal as long as the /3-coupling is 
finite. 

Analogously, by taking a scalar product of the rhs' of 
Eqs. (|l|) and (|l|) by ^^(t), it can be shown that 



in 



Wi tn- 



(28) 



Finally, we combine relations (pl|), (22), (p6|), and (12 
into a closed homogeneous equation with respect to f|: 
f| == U-^W-^UiWi f|. In order for this equation to have 
a non-trivial solution, the Jacobian matrix must satisfy 
the following condition: 



det{l - U^W^UiWi) =0. 



(29) 



The roots e^ of this equation determine the dispersion 



Ref. 



■R). 



Here, e^ 



-(0) 



as 



law of electrons (see, e.g 

L/la -^ 0. 

Alternatively, the relations ^, (JH), (H), and (|2 
could have been resolved with respect to g| , which would 
have led to the permuted Jacobian matrix UiWiU-^Wp 
However, such a permutation leaves the determinant un- 
changed and therefore the dispersion equation (ES) would 
have been still applicable. 

To obtain the solution of Eq. (E3) for given values of 
the external parameters luj/lai luj/lf3, and kyl^ we trun- 
cated the infinite matrix l — U-^W^UiWi down to the first 
N rows and columns and used a numerical root finder to 
obtain the zeros of its determinant. Owing to the con- 
veniently chosen bases (ISt) and (E9[), the roots of the 
equation (|29|) converge very quickly to their exact values 
as N is increased. For instance, for l^j/la, luj/lp ^ 1, 
it suffices to take iV « 30 to find the ten lowest energy 
levels with a net error < 10^^. 

The energy spectrum Sx = Sxikyl^i) calculated by this 
procedure is shown in Fig. 0. Let us first consider the case 
of neghgibly weak /3-coupling, when l^/lp -^ [Fig. ||(a)]. 
Here we see two-fold spin degeneracy of all quantum lev- 
els at fcj^ = 0. As soon as we move away from the point 



ky = 0, the SO interaction lifts the degeneracy and pro- 
duces an energy splitting A/j = ej, — e^ proportional to 
ky, where ej,-'- are energies associated with the "up" and 
"down" spin polarisations respectively. This i Juiya j ; b | 0- | 
haviour agrees with both theoretin^|pr|edictionsliilE3E3'C3 
and experimental observational3'EjO'E3 of the Rashba 
splitting in 2D systems. 
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FIG. 3. The transverse energy Sx vs. kylu> for finite 
Qf-coupling (L/L = 0.3): (a) L/'/3 = 0; (b) L/';3 — 0.1. 



The linear dependence of A^ on ky can be easily de- 
duced analytically from Eqs. (||) and (0) within first-order 
perturbation theory with respect to the parameter luj/la- 
Indeed, if we treat the rhs' of Eqs. (||) and (|^) as pertur- 
bations, then the first-order correction AeJ.J- = el^ — ei ' 
to the unperturbed energy (|l6|) will be given by 



Ael^ 



±xlP^'^P]^'^, 



A, 



2,/PpP™^ 



(30) 



(31) 



Here P"" 
operators P^i: 



^ I are diagonal matrix elements of perturbation 
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(32) 



The total (perturbed) wave function is given by $|| = 
(/)!|*i + J$? , . Fig. ^ shows graphs of <i>^| (t) for n = and 
n = 1. 



^" = i:(4 + ^-'" 



(33) 



To calculate the matrix elements P|j" we use the prop- 
erty (pG) and take into account that 



i'l^nl^W ='5™n 



for 



L/lfj = 



(34) 



[see Eqs. (H) - (||)]. As a resuh, we obtain Pjf ~ 
{luj/la){kylu) and hence the Rashba splitting A^ ( pi] ) is 
governed by 
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L/^/3 = 0, 



(35) 



e^°\n 



i.e. it is indeed proportional to the longitudinal wave 
number ky. 

The perturbative result ( p5| ) applies as long as A;? <C 
l)-ei°^(n), i.e. if (L/Za)(A:j,L) < 1. For 
^w/^Q = 0.3 this condition restricts kyl^ to values much 
less than about 3. In Fig. |3(a) we see that the linear 
behaviour of the energy splitting for all the curves holds 
true until kyl^ w 2 (dotted line). Within the region 
kyluj ^ 2 the asymptotic ( |35| ) gives a very good fit to all 
the energy branches. However, as soon as kyl^ becomes 
larger than fv 2.5, the dispersion curves start to bend. 
Eventually this leads to an anticrossing of branches cor- 
responding to quantum levels with neighbouring discrete 
numbers n. This fact drastically contrasts the simpler 
situation with no confining potential (ph, where the liar 
early split-spectrum (BS) represents the exact solutionc3. 

Usuallytll anticrossing is associated with similar sym- 
metry ( "hybridisation" ) of underlying states (wave func- 
tions). To reveal possible symmetries in our case, we 
return to the perturbation theory with respect to lui/la 
and consider first-order corrections S^^,{t) to the wave 
functions (EG): 



6^\{t)^Y. 
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(36) 



Owing to the orthogonality condition (p4) the second 
(constant) term in Eq. (|3^ ) does not contribute to S^^At) 
( pq ) for any m ^ n. At the same time, accordin g to 
Eqs. (|2(]| ) and (p4|), the first (differential) term in Eq. (|3^ ) 
gives rise to non-zero matrix elements Pmn for m = n+l 
and m = n— 1. As a result, S^^At) ( p6[) takes the form: 

for L/lp -^ 0. (37) 




FIG. 4. Unperturbed wave functions {l^/la ~ 0, solid 
curve) and wave functions modified by the a-coupfing 
{lui/la = 0.3, the dashed curve describes the "up" and the 
dotted curve the "down" spin orientations): (a) n = 0; (b) 
n=l. Zero /3-coupling {l^/lp — 0) is assumed. For definite- 
0. 



Eq. ( pTj ) demonstrates that the n-th quantum state is 
no longer independent of its "neighbours", i.e. of the 
(n ± l)-st states. In turn, the (n -I- l)-st state now de- 
pends on its own nearest neighbours, i.e. on both the 
n-th and the (n -I- 2)-nd states. Owing to this interstate 
coupling, all spinors <&" partially acquire the symmetry 
properties of spinors $"^1. Moreover, since the result 
( p7| ) is perturbative, it is plausible that the interstate 
coupling becomes more pronounced when the perturba- 
tive approach breaks down, i.e. when kyl^ ^ 2.5. This 
"hybridisation" of electron states accounts for the anti- 
crossing of neighbouring energy branches in Fig. @(a). 

After passing the relatively narrow anticrossing region 
2.5 < kylt^ < 4.5, all the curves in Fig. ||(a) straighten out 
and adopt a very accurate linear behaviour for kyluj ^ 5. 
There are two main features that attract attention in the 



infinite interval kyl,^ > 5: (i) there are no further an- 
ticrossings of approaching energy branches. Even those 
with closest quantum numbers n, which anticrossed at 
smaller values of kyl^^, now cross; (ii) all the straight lines 
el.^{n) go down (up) with the same slope « ^luj/la which 
is independent of a level number n. Although statement 
(i) means asymptotical smallness (rather than absolute 
absence) of anticrossings in the region of large kyli^ > 5 
in Fig. 0(a), these anticrossings are much weaker than 
any other anticrossings in the system and therefore can 
be left out. 

To explain these properties we should return to 
Eqs. (|), (0), and (^. According to Eq. ^), the spa- 
tial variation scale of functions </'ti(i) is of the order 
unity. This means that the derivatives $'|| on the rhs' 
of Eqs. (g) and (0) respectively can be roughly estimated 
as $11, at least for |t| ^ 1 and sufficiently small values 
of luj/la- Hence, it seems plausible that within the re- 
gion, where kyl^ ^ 1, the terms $'|| can be neglected 
in comparison with {kyl^)(^ ^-^ . Taking this into account 
together with the condition l^/lfs = 0, we reduce Eqs. (0) 
and (ItI) to 



$'/ 



n 



t^)^^{t)^{L/L){kyL)^i{t) (38) 



(e, - t2) $^(t) = {L/la.){kyL)^^{t). (39) 



Unlike the general problem (pi) and (M), Eqs. (BSh and 
( |39| ) can be simply decoupled by a unitary transforma- 
tion in the spin space. The corresponding matrix has the 
form 



l/\/2 -l/\/2\ 
1/V2 1/V2J 



(40) 



It is straightforward to verify that eigenvalues (energies) 
of the decoupled equations are 



^JJ = 2n+lTiL/la)ikyL) 



(41) 



and eigenfunctions are identical to those given by 
Eq. ( p^ at luj/lp = 0. The formula ( pl| ) yields the ex- 
pected linear dependence of both e], and e^ on kyl^j with 
slopes =f(Z^/Zq) independent of a level number n. The 
asymptotics (^Tj) give a very accurate fit to the spec- 
tral pattern in the region kyl^i ^ 5 in Fig. H(a). Since 
Eqs. ( p^ ) and ( |3^ ) can be decoupled by a simple rotation, 
states <&? and $? turn out to be intrinsically independent 
of each other as well as of all the other states. This fact 
alone explains why the anticrossing of energy branches is 
not observed for sufficiently large values of kyl^^. 

Let us now consider how the energy spectrum is modi- 
fied by switching on the J^-couphng [Fig. ||(b)]. From vi- 
sual comparison of Figs. y|(a) and (b) we see that the main 
effect of the /3-coupling is to enhance considerably the an- 
ticrossing of "neighbouring" spectrum branches. More- 
over, the strength of the anticrossing now depends on 
the quantum number n and grows as n is increased. An 



interesting consequence of this behaviour is an essential 
reduction of the linear Rashba energy splitting Ajj ex ky , 
contrasting the expectation that the additional mecha- 
nism of the SO interaction should intensify the splitting 
rather than suppress it. 

To understand t he pe culiarities of Fig. p|(b) we note 



that in subsection [I B it was shown that a finite /3- 
coupling leads to spatial separation between "up" and 
"down" spinor components (see Fig. |]). As a result, the 
orthogonality condition (|34|) no longer applies. Instead, 
the asymptotics (^21) - (|15|) should be used. If we now 



calculate the Rashba splitting An (pi|), then we will see 

|(^?r^ ) give non-zero (viz linear 



that scalar products ( (fuiyii 

in luj/lp) contributions to the diagonal matrix elements 
, so that /S.R is now described by [cf. Eq. (IsF 
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From this formula it immediately follows that the energy 
splitting /S.JI diminishes in comparison with Fig. 0(a) by 
an amount (n -t- l/2)(li^/lp) that is proportional to the 
quantum level number n. This conclusion agrees well 
with data presented in Fig. @(b). 

The result (B3) can be interpreted in the language of 
interstate coupling. Actually, it is easy to verify that 
the asymptotics (|l^ ) - (|5|) give rise to an additional 
(proportional to l^j/lfj) term ^^'^(i) in the wave function 
correction (5$','^ (i) (|^: 



S<^-{t)^l^f^f{kyL) 

Z La 1/3 



±xl'^{kyL)cf>'^ + Ht) 



+ I ^{n+l){n + 2)<t>]^+\t) 



i y;i(;r^0"r^(t) 



(43) 



This term, as well as (5$?i (p^) , involves "nearest" unper- 
turbed wave functions (/)?| . Since parameters luj/la and 
lu:/lf3 are independent of each other, we conclude that the 
/?-coupling enhances the hybridisation between an n-th 
and (n ± l)-st electron states and thereby the anticross- 
ing of corresponding energy branches. The strength of 
the hybridisation grows with the growth of n. 

Because of the enhanced interstate coupling, the an- 
ticrossing of neighbouring energy branches in Fig. g|(b) 
is not restricted to a narrow region of relatively small 
values of kyl,^ < 5 as it was in Fig. |3|(a). Instead, we 
can see that for all n > there exists a second anticross- 
ing in the region kyl^, ~ 7—10. Moreover, for larger 
quantum numbers n > 4 a third anticrossing emerges at 
kyL « 13 - 14. 



D. Ballistic conductance 

To apply the results of the previous subsections to the 
study of the effect of the SO interaction on the ballistic 
conductance of a long QIDES at low temperature, we 
must relate its conductance to its energy spectruni,-|Here 
we do this using the two-probe Landauer formulaca 



G 



Giep) = -j^M{eF). 



(44) 



where G is the ballistic conductance, ep is the Fermi 
energy, and Miep) is the number of occupied electron 
subbands which propagate in the same direction: 



^(^^) = E E E ^ [^^ - ^™"("' *)] 



(45) 



n i s=T,i 



Here e^j„(n, i) is the energy of the i-th minimum in the 
n-th electron subband with the spin orientation s. 9{x) 
is the Heaviside unit step function. As will be seen later 
the SO interaction can produce multiple minima in each 
ID subband owing to the anticrossings between different 
subbands. Since £^j„(n, i) can be found directly from 
the dispersion law of electrons, the conductance G (H) 
turns out to be completely defined by the energy spec- 
trum alone. Therefore, as long as the expression (H) 
applies, the knowledge of the electron energy levels in a 
QIDES will be sufficient to predict the behaviour of the 
ballistic conductance as a function of the Fermi energy. 

The applicability p£ Eq. (E3), as well as of the general 
scattering approachea to quantum transport, is essen- 
tially based on the condition that a current must travel 
in any ID electron subband without scattering into any 
other. In Refo it was shown that this current conserva- 
tion condition holds true in quite general circumstances 
including placing a QIDES in both a finite external mag- 
netic field and an arbitrary external electrostatic poten- 
tial. However, the proof in Refo ignored spin degrees of 
freedom and therefore did not take into account poten- 
tials acting in the spin space. Such neglect is not valid 
for our problem where the SO interaction Hamiltonian 
(0) leads to a highly non-trivial role for the spin in form- 
i ng th e energy spectrum of electrons (see subsections H A 

of Ref . 

conserved in the presence of a finite SO coupling. 

The first step is to define the matrix elements of the 
current density jmn (r) for the case where the wave func- 



j)_For that reason we cannot rely on the conclusion 
'.e3 but should check explicitly if the current is still 



tions ^„(„\(r) are spinors 
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(46) 



Here the dagger denotes the hermitian conjugate and 
p = —iKV. It is straightforward to verify that the di- 
vergence of the vector jmn(i') is given by 




VHI)^„-^1{vhA. (47) 



V •jm„(r) = 



We now suppose that the Hamiltonian of the system has 
the general form H = p'^/2M+Q, where Q is a hermitian 
operator (Q^ = Q) that acts in both coordinate and spin 
spaces. As applied to our problem, Q ~ Vlc + Hso [see 
Eq. (||)]. Using this Hamiltonian and the Schrodinger 
equation H^ ~ E"^, we express V^\I'J„ and V^^„ in 
terms of ^l^ and ^„ respectively and substitute the ex- 
pressions obtained into Eq. (47). As a result, we have 



V-j, 



,(r) = -(i?,„-i?„)*t„(r)vl'„(r) 



(48) 



where Em and En are energies corresponding to the 
states \l/m and 'i'n respectively. In the absence of any 
inelastic collisions, any scattering occurs between states 
of the same energy. So, without loss of generality, we 
can restrict ourselves to considering only equal energies 
Em = En in Eq. (Eq), in which case we find that 



V •j„m(r) = 



for 



E — E 



(49) 



This fundamental identity ensures local current conser- 
vation in the system. Once it has been established, if the 
system is translationally invariant in the longitudinal di- 
rection, further steps in proving the global conservation 
of current do not depend on specifics of rthe Hamiltonian 
and can be carried out in line with Ref.cZI. Therefore we 
arrive at the conclusion that the total current (i.e. cur- 
rent integrated over the cross section of the channel) be- 
tween states m and n at Em = En is equal to zero unless 
m = n. In other words, eigenstates of the Hamiltonian 
H are perfect current-carrying states that are free from 
scattering even in the presence of arbitrary SO coupling. 
For our problem this result implies that the spectrum 
of the Hamiltonian (0) is directly relevant to and com- 
pletely defines the ballistic conductance in the presence 
of the SO coupling. This allows us to use the simple Lan- 
dauer formula (E3), (|45|), in which the minima e^j„(n, i) 
of the energy subbands can be found from the analysis 
of the spectrum presented in subsections [I A - |H( 



Now we are in a position to discuss the features of 
the ballistic conductance in a QIDES subject to the SO 
interaction. For illustrative purposes we start with the 
"ideal" case of zero SO coupling. The corresponding sub- 
band energies e{n) = 2n + 1 + {kylu,Y are plotted in 
Fig. 11(a) (solid curves) as functions of kyl^. The de- 
pendence G{£f) can simply be deduced from this figure 
by moving a horizontal line e = ep from zero upwards 
and counting the number of points at which this line 
crosses the spectral parabolas. Since all the subbands in 
an ideal system are two-fold spin degenerate for any /cj,, 
this number coincides with the number M of propagating 
modes in the iQIDES. As a result, we restore the well- 
known picturcaU'H of ballistic conductance quantisation 



with equidistant jumps each of height 2e^//i [sohd curve 
in Fig. 1(a)]. 
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FIG. 5. The subband energy e in units of Tiijj/2 vs. kyli^: 
(a) zero SO coupling (solid curves), finite /3-coupling (dot- 
ted curves); (b) weak Q-coupling (luj/la < \/2); (c) strong 
Qf-coupling {luj/la > \/2) and zero /3-coupling; (d) a magnified 
bump on an n = energy branch in the anticrossing region; 
(e) strong a-coupling {luj/la > \/2) and finite /3-coupling. 
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FIG. 6. The conductance G in units of e^ /h vs. the Fermi 
energy £f in units of h,uj/2: (a) zero SO couphng and finite 
/3-coupling (solid line), weak a-coupling (dashed line); (b) 
strong Q-coupling and zero /3-coupling; (c) strong a-coupling 
and finite /3-coupling. 



Dotted curves in Fig. g(a) show subband energies 
£{n) = £x [n) + {kyluiY in the presence of the /3-couphng 
[see Eq. (^]. Each subband is two- fold spin degener- 
ate. It is seen that for any given value of ep the number 



M of forward propagating modes remains the same as 
it was in the ideal limit. Thus, we conclude that the (3- 
coupling alone does not affect the ballistic conductance. 
To detect the presence of the /3-coupling in a system one 
should exploit a unique feature that can only be due to 
the /3-coupling. As a reasonable example of such a fea- 
ture, we mention the spatial separation betwee n sp inor 



components that was established in subsection 11 B and 
illustrated in Fig. |l. Owing to this effect, the regions of 
the most probable distribution of electrons with opposite 
spin orientations will be shifted apart in the transverse 
direction. By analogy with ferromagnetic materials, this 
effect can manifest itself in the appearance of the mag- 
netisation of the electron system and its dependence on 
the transverse coordinate. 

Now we move on to the more complex case of fi- 
nite a-coupling. To start with let us imagine the fic- 
titious situation where the linear Rashba energy split- 
ting A;^ = 2{L/l^){kyL) [see Eq. (||) and Fig. |(a)] 
holds true not only for small values of kyl^ but for all 
kylc^. Here, the subband energies would be given by 
efl(n) = 2n -I- 1 -|- {kylui)^ ± {luj/la){kyl^). This formula 
describes parabolas of the same form as the solid curves 
in Fig. H(a) but shifted by an amount -ii{l/2){l^/la) along 
the abscissa and lowered vertically h\-iA / A){li^ / la)^ . The 
parabolas £||(n) were used in Ref.E3 to predict spon- 
taneous spin polarisation of ballistic electrons in quan- 
tum wires due to spin splitting. Whereas the real energy 
branches £f (n) and £i{n + \) anticross [see Fig. y(a)], the 
fictitious dispersion curves i-iin) and £|(n -I- 1) cross at 
the point kyl^, = (lu/la)^^- This point remains the same 
for all numbers n and goes to infinity as the SO coupling 
vanishes. It is easy to verify that for li^/la < V^ the 
crossing point kyl^ — {lui/la)^^ fies to the right of the 
point kyli^ = {l/2){la,/la) at which the parabola ex(n -1-1) 
has the minimum. At the same time, if l^j/la > v2, then 
the crossing point turns out to be to the left of the min- 
imum point. Below we will see that these simple conclu- 
sions following from the fictitious energy spectrum play 
an important role in determining the behaviour of the 
conductance G{eF)- 

Fig. ||(b) presents the energy subbands £]i{n) (solid 
and dashed curves) of a QIDES in the case of weak a- 
coupling when luj/la < V2 (more specifically, l^j/la = 
0.5). Dotted lines indicate the fictitious energy levels 
e|(n — 0) and £i{n = 1). In full accordance with the 
above conclusions, these levels cross to the right of the 
bottom of the parabola ej(n =1). As a result, the cross- 
ing angle turns out to be quite small. From Fig. ^(b) it 
is seen that this angle determines essentially the shape of 
the real energy curves e-\ {n = 0) and e^ (n = 1) within the 
anticrossing region 1-5 ^ kyl^i ^3- As long as the angle 
is small, the anticrossing remains very smooth and the 
curves Si{n = 0) and £i(n = 1) behave very much like 
the ideal parabolas in Fig. ^(a). The same observation 
is also true for higher quantum numbers n. Such simi- 
larity suggests that weak a-couphng (L/^q < \/2) does 
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not have a strong effect on the conductance. Indeed, it is 
easy to see by scanning Fig. ^(b) with the horizontal hne 
e = ep that the only effect of weak a-coupling on G{eF) 
is to shift the conductance quantisation steps down to 
lower energies by an amount {l/A){l^/la)'^ [see dashed 
lines in Fig. §|(a)]. 

As soon as the coupling constant lu)/la gets over the 
threshold of \/2, the spectral picture becomes much more 
interesting. Fig. |5|(c) shows a case when l^j/la — 1-8. 
Here the parabolas e'[{n — 0) and e|(n = 1) (dotted 
curves) cross to the left of the bottom of e j (n + 1 ) , which 
makes the crossing angle relatively large. The direct con- 
sequence of this is that the energy branches e^{n = 0) 
and e|(n = 1) cannot now anticross as smoothly as they 
did for lui/la < v2. Instead, in order to keep continu- 
ity, the lower curve [i.e. e^ {n = 0)] is forced to exhibit a 
non-monotonic portion ( "bump" ) within the anticrossing 
region. Although the height of the bump is rather small 
(~ 0.05 in the dimensionless units or "^A25 meV as de- 
termined from typical subband spacingala) , its geometry 
is a unique feature of the SO coupling mechanism [see a 
magnified plot of the bump in Fig. ra(d)] . The remarkable 
fact about it is that the bump contains a region where 
the electron energy decreases as the wavenumber grows. 
In other words, there exists an energy interval within 
which electrons are allowed to have a negative longitudi- 
nal group velocity Vy = h~^{de/dky) < for positive ky. 
Since this interval has a finite width and is surrounded 
by energy domains with positive group velocities, the 
curvature of the line e = e{ky) reverses sign within the 
negative- velocity interval and therefore there exists an in- 
flection point where d'^s/dky = [see Fig. |^(d)]. At this 

point the effective electron mass M = h [d^e/dkT) di- 
verges. Since the current passing through a ID subband 
is proportional to thn.product of the group velocity and 
the density of statesQij, the singular effective mass has no 
effect on the conductance in the absence of electron scat- 
tering. However, it could significantly change the depen- 
dence Giep) in the presence of a scattering mechanism 
(e.g., disorder or electron-electron interaction). 

As far as transport through ID electron subbands is 
concerned, it is clear from Fig. 0(d) that the negative 
propagating mode coexists with a forward propagating 
mode with the same spatial wave function. In this phys- 
ical situation it is likely that weak elastic scattering be- 
tweea these states would result in directed localisation 
Ref.EJ so that they would not be observed in conduc- 
tance. 

However, if both forward and backward electron modes 
contribute equally to the conductance, then this imme- 
diately brings us to the conclusion that the bump on 
the curve e|(n = 0) gives rise to a peak in the depen- 
dence G{eF) [see Fig. |g(b)]. The height of the peak is 
2e^//i and its width is defined by the height of the bump. 
Analogous bumps can be observed for all quantum levels 
n = 0, 1, 2, ... As a consequence, the sharp peak is seen 
in Fig. ra(b) not only on the first conductance plateau 



G — 2e^//i but also on the second one G — Ae^ /h and 
would still be seen on all further plateaus. 

Of course, the strong a-coupling {lui/la > V^), as well 
as the weak coupling (L/^a > V2), shifts the conduc- 
tance quantisation steps to lower values of the Fermi 
energy in comparison with the ideal situation where 
luj/la — 0. However, in contrast to the weak-coupling 
limit, strong coupling makes this shift much larger for 
the first step than for the others. This effect is clearly 
seen in Fig. p(b) as opposed to the ideal conductance 
quantisation shown by dotted lines. The larger shift of 
the first step is explained by the fact that the n — 
quantum state does not have a neighbour with the next 
lowest quantum number n. The energy level £i{n = 0) is 
therefore not forced to anticross with any other (lower ly- 
ing) energy branch and therefore nothing affects its linear 
(Rashba) dependence on ky [see Fig. ||(a)]. 

The final step of our analysis is to consider the case 
where both strong a-coupling {luj/la — 1-8) and rela- 
tively weak /3-coupling {luj/lp — 0.2) are present in the 
system. The corresponding energy bands are shown in 
Fig. 0(e). As we demonstrated in subsection |1IC| , fi- 
nite /3-coupling tends to suppress the energy splitting 
caused by the a-coupling [cf. Figs. ||(a) and ||(b) and 
see Eq. (p2|)]. Effectively, this suppression can be inter- 
preted as the enhancement of the anticrossing between 
neighbouring energy branches £|(«) and £i{n + 1). In 
Fig. ^(e) this effect is seen as a decrease in the gap be- 
tween energy branches etfn) and ex(n) {n — 1, 2, 3, . . .) 
in comparison with Fig. ^(c). Moreover, in accordance 
with Eq. (B3) , the gap becomes monotonically smaller as 
n grows. The enhanced anticrossing has a drastic effect 
on the energy bumps created by the strong a-coupling: 
all the bumps except for the lowest one are now smoothed 
away and do not give rise to the two additional (forward 
and backward) propagating modes. As a result, the con- 
ductance exhibits only one sharp peak that occurs on 
the plateau G = 2e^/h [see Fig. ^(c)]. A second effect of 
the enhanced anticrossing on the conductance is that the 
conductance quantisation steps (starting from the sec- 
ond one) are now located closer to the ideal steps than 
they were for zero /3-coupling. As each next step appears, 
its distance from the corresponding ideal step diminishes 
until eventually they fuse. 

There is another effect of the enhanced energy anti- 
crossing on the conductance that is not seen in Figs. ^ 
but could possibly be detected in. the .presence of impu- 
rities. It is found experimentallyEilE3 that the ballistic 
conductance quantisation breaks down as the constric- 
tion becomes too long (at least > 2 /im). This effect 
can be attributed to potential fiuctuations caused by the 
rauriom distributions of remote impurities. Zagoskin et 
aO showed analytically that the degradation of quan- 
tised conductance decreases exponentially with the ratio 
of the ID subband spacings to the standard deviation of 
impurity potential fluctuations. Obviously, the subband 
spacings depend on the presence of the SO interaction 



12 



(see Figs. ^. By comparing Figs. ||(c) and (e) we see that 
a finite /3-coupling pushes neighbouring energy branches 
apart and thus increases the subband spacings (possibly 
by tens of per cent) . This means that the presence of the 
/3-mechanism of the SO interaction could allow the ex- 
perimental observation of the conductance quantisation 
structure in longer samples. 



III. CONCLUSION 

We have studied theoretically the influence of the 
spin-orbit (SO) interaction on the energy spectrum and 
the ballistic conductance of a quasi-lD electron system 
(QIDES) formed by lateral electric confinement of a 2D 
electron gas. The presence of the confining potential 
proves to be crucial in two ways. First, it alters sig- 
nificantly the effect of the quantum well asymmetry on 
the energy spectrum in comparison w ith a purely 2D sys- 
tem. As it was shown by Rashba et acS'Ea, in 2D electron 
gases this asymmetry gives rise to SO coupling (which 
we refer to as the a-coupling) that manifests itself in a 
lifting of the spin degeneracy of electronic states. The 
accompanying energy splitting Ar was found to be pro- 
portional to the in-plane electron wave number k. We 
have demonstrated that in a QIDES the function Afl(fc) 
is non-monoinnic [Fig. ^(a)] with the standard linear de- 
pendenceE3'E3 Ajj ex k being observed for relatively small 
and large values of k only. Such a drastic change in 
behaviour is explained by the essentially different effect 
of SO coupling operator P|| (33) on unperturbed wave 
functions depending on whether or not the lateral confin- 
ing potential is taken into account. Indeed, in the zero- 
potential case the unperturbed wave functions are plane 
waves characterised by the continuous in-plane wave vec- 
tor k and the spin. The action of the operator P|| on 
any such wave function reduces to a simple renormalisa- 
tion of its amplitude with the renormalising factor be- 
ing independent of coordinates. As a consequence of this 
simple effect, wave functions belonging to different values 
of k remain independent of each other and the SO cou- 
pling manifests itself in linear (monotonic) energy split- 
ting Afl; c>c k which is the same for all quantum states. 

In contrast to this, the wave functions in the presence 
of a finite confining potential have a more complicated 
structure (HG) and are characterised by the discrete quan- 
tum number n instead of the continuous transverse wave 
number. According to Eq. (pOl), the action of the oper- 
ator P|| on the ri-th wave function includes projecting 
it onto states with the next (preceding) closest quantum 
numbers n±l. In turn, the (n±l)-st states are projected 
onto their "neighbours" with numbers n and n±2 respec- 
tively. As a result, the operator P-^i couples effectively 
any n-th and (n ± l)-st wave functions. In other words, 
once the SO coupling has been taken into account, the 
n-th and {n ± l)-st wave functions cease to be indepen- 
dent and possess symmetry elements of each other. This 



partial symmetry between states leads to an anticrossing 
of the closest (neighbouring) energy branches in Fig. |^(a) 
and hence to the non-monotonic dependence Afl(fc). 

Apart from the interplay with the familiar (quantum- 
well-asymmetry or Rashba) mechanism of the SO inter- 
action, the lateral confining potential by itself appears to 
be a source of additional dynamical coupling between the 
orbital and spin degrees of freedom of an electron. This 
coupling (which we refer to as /3-coupling) originates from 
the natural spatial non-uniformness of the confining po- 
tential. A typical variation scale of this potential lies 
within the wide range '^ 10 - 1000 nm, which makes the 
accompanying electric field sufficiently strong to compete 
with the quantum-well-asymmetry field. This competi- 
tion may become especially noticeable in square quan- 
tum wells with a relatively weak Rashba contribution. 
Whereas the quantum-well field is normal to the device 
plane, the confinement-induced electric field is parallel 
to the plane. It is also spatially non-uniform as long as 
the coordinate dependence of the confining potential is 
more complex than linear. These features make the /3- 
coupling an essentially different (new) mechanism of the 
SO interaction in a QIDES which cannot be taken into 
account by simply adjusting the Rashba interaction con- 
stant. This claim is confirmed by Fig. 0(b) that demon- 
strates the combined effect of both the Rashba and the 
confinement-induced couplings on the energy spectrum. 
From comparison of Fig. ||(b) with Fig. y(a) we see that 
the major role of /3-coupling is to reduce the Rashba en- 
ergy splitting Aji (see region kyl^ < 2) rather than to 
give a positive correction to it. It is crucial that the 
reduction depends on the quantum number n and mono- 
tonically grows as n increases. This fact is a clear man- 
ifestation of the position dependence of the electric field 
created by the confining potential and indicates the in- 
dependent nature of /3-coupling. The suppression of the 
Rashba energy splitting is a part of the overall effect of 
the /3-coupling on the underlying electron system which 
consists of a considerable enhancement of the hybridis- 
ation between different quantum states and eventually 
leads to a more pronounced anticrossing of neighbouring 
energy branches. 

In addition to the investigation of the SO-effects in 
the energy spectrum of electrons, we have discussed the 
possible manifestations of the SO interaction in the bal- 
listic conductance of a QIDES. The key point of our ap- 
proach to the conductance is the fundamental current- 
conservation identity that was proven in subsection II D . 



According to this identity, the electron eigenstates that 
were found as the solution to the spectral problem are 
perfect current-carrying states. A current can travel in 
any of these states without scattering into any other. 
This property therefore allows the ballistic conductance 
to be calculated directly from the energy spectrum with 
the help of simple Landauer formula (H). 

An analysis of the ballistic conductance G reveals that 
the /3-coupling alone does not affect the dependence of 
G on the Fermi energy ep- This fact is illustrated by 
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Fig. ^(a) where we see that the /3-coupling reduces the 
curvature of the parabohc energy bands (cf. dotted 
and sohd curves) , while the band edges (bottoms of the 
parabolas) remain anchored. Thus, /3-coupling cannot be 
detected by measuring the ballistic conductance in sys- 
tems where the SO interaction appears only in a form of 
the /3-coupling (e.g. in square wells). Instead, one should 
use experimental methods that allow direct observation 
of the electron energy spectmjn in a QIDES, e.g. mag- 
netotunneling measurementsCj. 

In the presence of a-coupling the behaviour of G{ep) 
is essentially determined by the strength of the SO inter- 
action. If the a-coupling is not too strong [luj/la < V^), 
then its only effect on the conductance will be shifts of 
the conductance quantisation steps to lower Fermi en- 
ergies in comparison with an ideal (i.e. with zero SO 
interaction) situation. We note that such shifts are not 
related to the lateral confining potential and should also 
be present in purely 2D systems. This effect should be 
deferable in transconductance measurements (see, e.g., 
Ref.El) which determine both the conductance and the 
subband spacings simultaneously. 

In the limit of strong a-coupling {luj/la > V^) the 
electron energy bands take on a very interesting form 
[Fig. Mc)]. The most remarkable feature is the appear- 
ance of narrow energy intervals where two additional (for- 
ward and backward propagating) electron modes exist 
[Fig. [^(d)]. Such intervals can be found in each ID sub- 
band, starting from the lowest one. The additional elec- 
tron modes have similar magnitude group velocities but 
propagate in opposite directions. They have almost iden- 
tical subband wave functions and therefore would be sus- 
ceptible to strong intermode scattering in the presence of 
disorder. However, in a sufficiently pure QIDES, the ad- 
ditional electron modes give rise to the unusual periodic 
sharp steps in G{eF) shown in Fig. ||(b). This picture is 
changed by switching on relatively weak {1^/1/3 <^ l^j/la) 
/3-coupling. As we mentioned above in the discussion of 
the energy spectrum, the /3-coupling enhances the anti- 
crossing of energy levels initiated by the a-coupling. As 
applied to the conductance, this effect leads to quench- 
ing the sharp conductance peaks by the /3-coupling. The 
existence of the single peak (or just a few of peaks) in 
the dependence of G on ep could be a clear experimen- 
tal indication of the presence of the /3-coupling in the 
system. 

The most crucial point for the experimental observa- 
tion of the conductance peaks in Fig. J6| is to make the 
ratio luj/la sufficiently large {lu/la > V^- In typical sys- 
tems, where enhancing the SO interaction is not paid spe- 
cial attention, the value of lui/la hardly exceeds 0.5. An 
additional, at least threefold, increase in lui/la could be 
achieved by using: (i) materials with light carrier masses; 
(ii) strong (narrow) lateral potential confinement ^ 100 
nm; (iii) heterojunctions (triangular well) rather than 
quantum- well heterostructures (square well); (iv) a back 
gate voltage to maximise the interface (Rashba) electric 
field. 



We hope that the results presented in this paper will 
stimulate further experimental and theoretical work with 
the aim of understanding the role of the spin-orbit inter- 
action in determining the transport properties of quasi- 
ID systems. 



ACKNOWLEDGMENTS 

AVM thanks the ORS, Cambridge Overseas Trust and 
Corpus Christi College for the Research Studentship. 
CHWB thanks the EPSRC for the Advanced Fellowship. 
This work was funded by the EPSRC. 



^ T. J. Thornton, M. Pepper, H. Ahmed, D. Andrews, 

G. J. Davies, Phys. Rev. Lett. 56, 1198 (1986). 
^ C. W. J. Beenakker and H. van Houten Quantum trans- 
port in Semiconductor Nano structures, vol. 44 of Sohd state 
physics (Academic Press, New York, 1991). 
M. J. Kelly Low- dimensional semiconductors: material, 
physics, technology, devices (Oxford University Press, Ox- 
ford, 1995). 

B. J. van Wees, H. van Houten, C. W. J. Beenakker, 
J. G. Williamson, L. P. Kouwenhoven, D. van der Marel, 
and C. T. Foxon, Phys. Rev. Lett. 60, 848 (1988). 
^ D. A. Wharam, T. J. Thornton, R. Newbury, M. Pepper, 
H. Ahmed, J. E. F. Frost, D. G. Hasko, D. C. Peacock, 
D. A. Ritchie, and G. A. C. Jones, J. Phys. C 21, L209 
(1988). 
^ K. J. Thomas, J. T. NichoUs, M. Y. Simmons, M. Pepper, 
D. R. Mace, and D. A. Ritchie, Phys. Rev. Lett. 77, 135 
(1996). 
"^ A. Kristensen, J. Bo Jensen, M. Zaffalon, C. B. S0rensen, 
S. M. Reimann, M. Michel, A. Forchel, J. Appl. Phys. 83, 
607 (1998). 
** K. J. Thomas, J. T. NichoUs, N. J. Appleyard, M. Y. Sim- 
mons, M. Pepper, D. R. Mace, W. R. Tribe, and 
D. A. Ritchie, Phys. Rev. B 58, 4846 (1998). 
^ S. Datta Electronic transport in mesoscopic systems (Cam- 
bridge University Press, Cambridge, 1995). 

^'^ V. K. Thankappan Quantum Mechanics (John Wiley & 
Sons, New York, 1993); L. I. Schiff Quantum Mechanics 
(McGraw-Hill, New York, 1968). 

^^ L. D. Landau, E. M. Lifshitz Quantum Mechanics (Perga- 
mon Press, Oxford, 1991). 

^2 G P. Fisher, Am. J. Phys. 39, 1528 (1971). 

^^ H. L. Stormer, Z. Schlesinger, A. Chang, D. C. Tsui, 
A. C. Gossard, W. Wiegmann, Phys. Rev. Lett. 51, 126 
(1983). 

^■^ D. Stein, K. von Klitzing, G Weimann, Phys. Rev. Lett. 
51, 130 (1983). 

^^ B. Das, D. C. Miller, S. Datta, R. Reifenberger, 
W. P. Hong, P. K. Bhattacharya, J. Singh, and M. Jaffe, 
Phys. Rev. B 39, 1411 (1989). 



14 



^^ B. Das, S. Datta, R. Reifenberger, Phys. Rev. B 41, 8278 

(1990). 
^^ P. D. Dresselhaus, C. M. A. Papavassiliou, R. G. Wheeler, 

and R. N. Sacks, Phys. Rev. Lett. 68, 106 (1992). 
^* B. Jusserand, D. Richards, H. Peric, and B. Etienne, Phys. 

Rev. Lett. 69, 848 (1992). 
^^ G. Goldoni and A. Fasohno, Phys. Rev. Lett. 69, 2567 

(1992). 
'^° G. L. Chen, J. Han, T. T. Huang, S. Datta, D. B. James, 

Phys. Rev. B 47, 4084 (1993). 
^^ J. Nitta, T. Akazaki, H. Takayanagi, Phys. Rev. Lett. 78, 

1335 (1997). 
^^ G. Dresselhaus, Phys. Rev. 100, 580 (1955). 
^^ J. Luo, H. Munekata, F. F. Fang, P. J. Stiles, Phys. Rev. 

B 41, 7685 (1990). 
'^'^ R. Eppenga, M. F. H. Schuurmans, Phys. Rev. B 37, 10923 

(1988). 
^'^ A. G. Mal'shukov, K. A. Chao, M. Willander, Phys. Rev. 

Lett. 76, 3794 (1996). 
2^ V. M. Edelstein, J. Phys.: Condens. Matter 7, 1 (1995). 
^■^ T. Ando, A. B. Fowler, F. Stern, Rev. Mod. Phys. 54, 437 

(1982). 
^** E. L Rashba, Fiz. Tverd. Tela (Leningrad) 2, 1224 (1960) 

[Sov. Phys. - Solid State 2, 1109 (I960)]. 
^^ Yu. A. Bychkov, E. I. Rashba, Pis'ma Zh. Eksp. Teor. Fiz. 

39, 66 (1984) [JETP Lett. 39, 78 (1984)]. 

^° S. I. Dorozhkin, E. B. Ol'shanetskii, Pis'ma Zh. Eksp. Teor. 

Fiz. 46, 399 (1987) [JETP Lett. 46, 502 (1987)]. 
^^ T. Hassenkam, S. Pedersen, K. Baklanov, A. Kristensen, 

C. B. Sorensen, P. E. Lindelof, F. G. Pikus, G. E. Pikus, 

Phys. Rev. B 55, 9298 (1997). 
^2 A. Yacoby, H. L. Stormer, K. W. Baldwin, L. N. Pfeiffer, 

and K. W. West, Solid State Commun. 101, 77 (1997); 

A. Yacoby, H. L. Stormer, N. S. Wingreen, L. N. Pfeiffer, 

K. W. Baldwin, and K. W. West, Phys. Rev. Lett. 77, 4612 

(1996). 
^^ S. Datta and B. Das, Appl. Phys. Lett. 56, 665 (1990). 
^* B. E. Kane, Nature 393, May 14 133 (1998). 
^^ G. Fasol and H. Sakaki, Phys. Rev. Lett. 70, 3643 (1993); 

Appl. Phys. Lett. 62, 2230 (1993). 
'^^ V. V. Kveder, Yu. A. Osip'yan, and A. I. Shalynin, 

Pis'ma Zh. Eksp. Teor. Fiz. 40, 10 (1984) [JETP Lett. 

40, 10 (1984)]; V. V. Kveder, V. Ya. Kravchenko, 



T. R. Mchedlidze, Yu. A. Osip'yan, D. E. Khmel'nitskii, 
and A. I. Shalynin, Pis'ma Zh. Eksp. Teor. Fiz. 43, 202 
(1986) [JETP Lett. 43, 25 (1986)]. 

^^E. I. Rashba, Usp. Fiz. Nauk 84, 557 (1964) [Sov. Phys. 
Usp. 7, 823 (1965)]. 

^* S. E. Laux, D. J. Frank, F. Stern, Surf. Sci. 196, 101 (1988). 

^^ H. Drexler, W. Hansen, S. Manus, J. P. Kotthaus, M. Hol- 
land, S. P. Beaumont, Phys. Rev. B 49, 14074 (1994). 

*''B. Kardynal, C. H. W. Barnes, E. H. Linfield, 
D. A. Ritchie, J. T. NichoUs, K. M. Brown, G. A. C. Jones, 
and M. Pepper, Phys. Rev B 55, R1966 (199 7). 

■*^ G.-H. Chen, M. E. Raikh, |cond-mat/9812085 . 

B. K. Ridley Quantum processes in semiconductors 
(Clarendon Press, Oxford, 1993). 

*^ G. E. Marques, L. J. Sham, Surf. Sci. 113, 131 (1982). 

** M. Abramowitz, I. A. Stegun Handbook of mathematical 
functions (Dover Publications, New York, 1965). 

*^Yu. A. Bychkov, V. I. Mel'nikov, E. I. Rashba, 
Zh. Eksp. Teor. Fiz. 98, 717 (1990) [Sov. Phys. JETP 71, 
401 (1990)]. 

^^ R. Landauer, IBM Journ. Res. Development 1, 223 (1957); 
R. Landauer, Philos. Mag. 21, 863 (1970); R. Landauer, in 
Localization, Interaction, and Transport Phenomena edited 
by G. Bergmann and Y. Bruynseraede (Springer- Verlag, 
New York, 1985). 

"■^ H. U. Baranger, A. D. Stone, Phys. Rev. B 40, 8169 (1989). 

** S. J. Koester, C. R. Bolognesi, M. Thomas, E. L. Hu, 
H. Kroemer, M. J. Rooks, Phys. Rev. B 50, 5710 (1994). 

''^ C. Barnes, B. L. Johnson, and G. Kirczenow, Phys. Rev. 
Lett. 70, 1159 (1993). 

^° K. J. Thomas, M. Y. Simmons, D. R. Mace, M. Pepper, 
and D. A. Ritchie, Appl. Phys. Lett. 67, 109 (1995). 

^^ G. Timp, R. Behringer, S. Sampere, J. E. Cunningham, 
and R. E. Howard, in Proceedings of the International Sym- 
posium on Nanostructure Physics and Fabrication, edited 
by M. A. Reed and W. P. Kirk (Academic, New York, 
1989), p. 331. 

^^ S. J. Koester, B. Brar, C. R. Bolognesi, E. J. Caine, A. Pat- 
lach, E. L. Hu, H. Kroemer, M. J. Rooks, Phys. Rev. B 53, 
13063 (1996). 

^^ A. M. Zagoskin, S. N. Rashkeev, R. I. Shekhter, and 
G. Wendin, J. Phys. Condens. Matter 7, 6253 (1995). 



15 



